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AN ANALYTICAL ANALYSIS OF THE DISPERSION PREDICTIONS 
FOR EFFLUENTS FROM THE SATURN V AND SCOUT -ALGOL 1 1 1 

ROCKET EXHAUSTS 


J. Briscoe Stephens 
Michael Susko 
John W. Kaufman 
C. Kelly Hill 


SUMMARY 


The NASA/MSFC Multilayer Diffusion Model is used to predict 
the dispersion of effluents from the rocket motors during the launch 
of Apollo 16 on April 16, 1972 at Kennedy Space Center and for a Scout- 
Algol III vehicle launched on August 13, 1972 at Wallops Island, 

Virginia. Predictions for two launches were made in real time to gain 
experience in making such predictions under operational conditions. The 
analyses of the dispersion of the rocket exhaust effluents from the 
Saturn V vehicle provide concentration isopleths for large liquid rocket 
motors. On the other hand, analyses of the effluents resulting from 
the launch of the Scout -Algol III show ground level concentrations from 
a small solid rocket motor. 

The analyses indicated that ground level concentrations of potential 
toxic effluents resulting from both launches are well below existing 
primary and secondary ambient air quality standards. 


DEFINITION OF SYMBOLS & TERMS 


Symbol 


Definition 

F 

= 

buoyancy flux 

C 

= 

if \ 

F 

= 

initial buoyancy term 

I 


3gQ 



I 

4c TnP 
P 

FM 

= 

percentage by 'weight of pollutant material in 
the fuel from Table 1 

H ,H 

= 

respective heat contents of liquid and solid fuels 

L S 

L 

= 

depth of the surface mixing layer 

M 

= 

molecular weight 

P 

= 

ambient pressure (mb) 

P{ z } 
1 TK J 


integral of the Gaussian probability function 
between minus infinity and the top of the K-th 
layer z 

TK 



P hrH 

Q 

= 

total weight of exhaust products in the stabilized 
exhaust cloud 


= 

( Qr )C r { z *i}) M 

Q 



rate of heat released by burning fuel 

F 


= 

H o W + H o W 
L L S S 


2 


Symbol 

Q 

K 

q r 


w ,w 

L S 


z 

c 

f 

g 

r 


R 

s si 

t a’z 


t* 


Definition 

effective heat released (cal) 

source strength in units of mass per unit depth 
of the K-th layer 

fuel expenditure rate from Table 1 
ambient air temperature (°K) 

respective fuel expenditure rates for liquid 
and solid fuel 

height above ground of any selected layer 
specific heat of air at constant pressure cal/(°Kg) 

2 

gravitationa 1 acceleration (9.8 m/sec ) 
initial cloud radius at the surface 

stability parameter 
time of layer breakdown 


H 


time required for the cloud to reach the 
stabilization height 


t z 
R ml 


u 


time after ignition in seconds 

time in seconds required for the vehicle to reach 

the height z of maximum rise of the ground cloud 
ml 

(obtained from Equation 1) 
mean wind speed 


z 


height of stabilized cloud 


z 


t 


midpoint of the K-th layer 

fz + z \/l 
\ BK TK }/ 
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Symbol 


Definition 


z = height of the base of the K-th layer 

BK 


z = height of the base of the L-th layer 

BL 


z 

L 


height in the L-th layer at which the concentration 
is calculated 


z 

mC 


maximum height of cloud rise for a continuous 
source 


z = maximum rise for an instantaneous source 

ml 


z = height of the top of the K-th layer 

TK 


z = height of the top of the L-th layer 

TL 


z = altitude above the pad in meters 

R 


y = entrainment constant 

C 

y - entrainment constant 

I 


ct = standard deviation of the concentration 

distribution of the stabilized ground cloud 

standard deviation of the alongwind concentration 
distribution in the K-th layer at distance x 


a 

xK 


a 

xLK 


a 


{K} 


xo 


CJ 

yo 


Ik] 


standard deviation of the alongwind concentration 
distribution in the L-th layer for the source 
originating in the K-th layer 

standard deviation of the alongwind concentration 
distribution in the K-th layer at cloud 
stabilization 

standard deviation of the crosswind concentration 
distribution in the K-th layer at cloud s tabilization 
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Symbol 

CT {K} 

zo 


a 


yK 


a 

yLK 


a zLK 


P 

59 

5z 


X 

P 


Centerline : 


Concentration ; 


Dosage : 


Ground Cloud : 


Plume Cloud : 


Definition 


standard deviation of the vertical concentration 
distribution in the K-th layer at cloud 
stabilization 

standard deviation of the crosswind concentration 
distribution in the K-th layer at distance x 

standard deviation of the crosswind concentration 
distribution in the L-th layer for the source 
originating in the K-th layer 

standard deviation of the vertical concentration 
distribution in the L-th layer for the source 
originating in the K-th layer 

density of ambient air (g/m ) 

vertical gradient of ambient potential tempera- 
ture 

peak or centerline concentration 


TERMS 


The radial vector in the direction of the mean 
wind direction whose origin is the launch site. 

is the amount of the effluent present at a specific 
time. The average concentration is the average 
amount present during the event. 

is the measure of the total amount of effluent 
(time integrated concentration) due to the 
vehicle launch at a specific location. 

That cloud of rocket effluents emitted during 
the initial phase of vehicle launch. This cloud 
has an ellipsoid shape. 

The cloud of rocket effluents emitted from th 
vehicle in flight. This cloud has a cylindrical 
shape whose height is defined by the vertical 
thickness of the layer. 
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TERMS 


Potential 
Temperature (e) 

Quasiadiabatic 

Layer : 

Stable Layer : 


is the temperature a volume of dry air would 
have if brought adiabatically from its initial 
state to the standard pressure of lOOOmb . 


is a layer in which the vertical potential 
temperature gradient is zero or less. 

is a layer in which the vertical potential 
temperature gradient is positive. 
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SECTION I. INTRODUCTION 


The NASA /MS FC Multilayer Diffusion Program has been utilized 
to predict the dispersion of the effluents from the booster during the 
launches of the Apollo 16 (Apollo Saturn 511) on April 16, 1972 at 
Kennedy Space Center and of a Scout vehicle launched on August 13, 1972 
at Wallops Island, Virginia. 

The National Environmental Policy Act of 1969 and the April 23, 
1972, guidelines of the Council on Environmental Quality require 
statements covering proposed federal actions that might affect the 
environment. Such statements are required for assessing the environ- 
mental impact of the Space Shuttle and other NASA space vehicle rocket 
motor effluents. Development of quantitative procedures for estimating 
the space vehicle rocket motor toxic fuel hazard has been underway for 
over a decade at the NASA Marshall Space Flight Center. Computerized 
and graphical procedures for estimating toxic fuel hazards have been 
developed [l] utilizing a Lagrangian model for estimations in the 
troposphere. In addition to estimates of atmospheric transport, dis- 
persal and decay of all airborne toxic material released as a result 
of normal launch operations, estimates must also be provided for cases 
involving fuel spillage, vehicle abort, or vehicle destruct situations. 

Universally accepted and adequately validated prediction 
techniques for the rocket motor effluent problem are not available, 
and much uncertainty exists concerning very important aspects of the 
problem, such as the amount and composition of the rocket engine 
effluents, and their dispersal and transport in the atmosphere. The 
available atmospheric measurements to ascertain the reliability of the 
description of rocket effluent dispersion models in the atmosphere are 
sparse and of questionable accuracy. On the other hand, the require- 
ments for estimating toxic fuel hazards clearly exist in order to 
establish special constraints on operations, test, and launch activities 
to assure allowable concentrations of these effluents will not be 
exceeded. The need for implementing a program for monitoring rocket 
engine exhaust effluents has been recognized for many years. As a 
result of informal discussions between representatives of NASA Head- 
quarters, Marshall Space Flight Center, Langley Research Center, and 
Kennedy Space Center, it became apparent that a NASA inhouse rocket 
engine effluent prediction and measurement program was desirable, 
possible, and feasible. 

The announcement by the President on January 5, 1972, that the 
United States should proceed at once with the development of the Space 
Shuttle Program resulted in renewed interest in the dispersion of 
potentially toxic engine effluents. The subsequent decision by NASA 


7 


to power the first stage of the Shuttle by Solid Rocket Motors added 
significant importance to the problem, since their major exhaust 
constituents contain large amounts of carbon monoxide (CO), hydrogen 
chloride (HC1) , and aluminum oxide (Al^O^) , all of which may be classi- 
fied as air pollutants and toxic. 

To initiate the program of predicting and measuring downwind 
concentrations of exhaust effluents, such predictions and measurements 
of the ground cloud and upper altitude plume generated by the Saturn V 
vehicle and, subsequently, a Scout vehicle were made. Our objective 
is to provide lead information for Solid Rocket Motors Programs whereby 
the various investigators may gain an insight into the behavior of the 
transport of the toxic exhaust effluents, and to improve the accuracy 
of various empirical parameters utilized in our dispersion calculations. 
To this end, we have selected two different actual examples — a large 
liquid rocket engine, the Saturn V, and a small solid rocket motor, the 
Scout -- to provide an initial baseline for the investigations of the 
dispersive effluent transport. 

The various descriptions for the models of dispersive transport 
of the rocket exhaust available with the NASA/MSFC Multilayer Diffusion 
Program required, because of the various phases of plume development 
and changing environment conditions, are summarized in the next section. 
The detailed formulations for the algorithms employed in the models 
along with a toxicity criteria are afforded in the appendix. Section 
III and IV give a detailed summary of the effluent transport predictions 
for the launch of Apollo 16 and for the launch of a Scout vehicle, 
respectively. 
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SECTION II. NASA/MSFC MULTILAYER DIFFUSION MODEL 


The spatial description, in terms of concentration and dosage, 
of the dispersive transport of effluents from a discrete source is 
afforded by the NASA/MSFC Multilayer Diffusion Model* Specifically, 
this application of the model is for the prediction of toxicity distri- 
bution associated with the rocket exhaust effluents emitted during the 
launch of a space vehicle in order to assess the resulting environmental 
impact. The dispersive description accorded by the Multilayer Diffusion 
Model is initiated at the point where the cloud of effluents reaches 
thermodynamic equilibrium with the environment and therefore depends 
strongly on the kinematic and thermodynamic profiles of the atmospheric 
conditions along with a knowledge of the exhaust effluents present in 
the cloud. 


The initial considerations in this section are given to the 
techniques of establishing the spatial location of the ground cloud 
equilibrium (see definitions). Secondly, a general discussion of the 
Lagrangian dispersion of a point source is given. The final discussion 
in this section explains how the Multilayer Diffusion Model incorporates 
the general diffusion description to account for the two stages of 
exhaust emission and accounts for the environmental effects. The 
significant mathematical expressions supporting these discussions have 
been included in Appendix I, II, and III. 

a. Altitude of Cloud Equilibrium 

The effluent cloud rise relations are employed to deter- 
mine at what altitude the ground cloud reaches equilibrium with the 
environment. The importance of this location is that it serves as the 
origin of the dispersive description. This equilibrium point is chosen 
as the origin in order to eliminate complex thermodynamic considerations 
and to limit the problem to solely kinematics. 

The burning of rocket engines results in the formation of 
a cloud of hot exhaust products which subsequently rises and entrains 
ambient air until an equilibrium with ambient conditions is reached. 

For normal launches, this cloud is formed principally by the forced 
ascent of hot turbulent exhaust products that have been deflected 
laterally and vertically by the launch pad hardware and the ground 
surface. The height at which this ground cloud stabilizes is determined 
by the vehicle type and atmospheric stability. The vehicle type deter- 
mines whether a continuous or instantaneous source model is required. In 
the instantaneous source model, spherical entrainment is assumed; that 
is, t\ie entrained ambient air enters the exhaust cloud uniformly from 
all directions. In the continuous source model, cylindrical entrainment 
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is assumed; that is, the entrained ambient air enters the cloud 
uniformly only on the sides of the cylinder and not the ends. Thus, 
this terminology -- continuous or instantaneous source — in reference 
to the cloud rise model does not imply the duration of the exhaust 
cloud, as it does in the diffusion model, but only implies the form 
of the entrainment process. The entrainment process is a function of 
the residence time of the vehicle on the pad. Experience to date 
indicates that the buoyant rise of exhaust clouds from normal launches 
of solid-fue led and small liquid-fueled vehicles is best predicted by 
using a cloud rise model for instantaneous sources; the cloud rise for 
large liquid -fueled vehicles is best predicted by the use of a cloud 
rise model for continuous sources . While no cloud rise data are 
available for on-pad aborts, cloud rise data from static tests of 
liquid-fueled rockets indicates that the use of a cloud rise model for 
continuous sources is appropriate. 

Each of the models forvcloud height is subdivided into 
two categories to account for the atmospheric temperature lapse rate. 

The model assumes that the atmosphere is either quasi-adiabatic or 
stable. Here the quasi-adiabatic is where the adiabatic atmosphere is 
the limit, which means that the potential temperature difference (A0) 
is zero or less, where the potential temperature difference is given by 

" 0 max cloud height" 9 surf ace • If this potential temperature differ- 
ence is positive, then the atmosphere is treated as stable. Since in 
most cases of interest there will be an inversion layer present, the 
stable cloud rise formula is the normally utilized relation. 

b. Generalized Diffusion Model 

The generalized diffusion model describes the kinematic 
transport — in terms of the temporal and spatial levels of concen- 
tration and dosage -- of the exhaust effluent constituents assuming 
the effluent cloud is in equilibrium with the environment. A Lagrangian 
model is assumed, where volumetric cloud expansion is about a reference 
point in a homogeneous fluid. For diagnostic and interpretation 
flexibility, this model is formated in a modular form for both concen- 
tration and dosage. 

The generalized concentration model for a nearly instan- 
taneous source is expressed as the product of five modular terms: 

Concentration = {Peak Concentration Term} x [Alongwind Term} x 
[Lateral Term} x [Vertical Term} x 
{Depletion Term} ; 
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whereas, the generalized dosage model for a nearly instantaneous source 
is defined by Lae product of four modular terms: 

Dosage = [Peak Dosage Terms} x [Lateral Term} 
x [Vertical Term} x [Depletion Term} 

Thus, the mathematical description for the concentration 
and dosage models permit flexibility in application to various sources 
and for changing atmospheric parameters while always maintaining a 
rigorous mass balance. 


Two obvious differences exist. First, the peak concen- 
tration term refers to the concentration at the point x, y = 0, z = H 
(where x is the wind direction and H is any height) and is defined by 
the expression 


/O v 3 / 2 

( 2 tt) a a a 
x y z 

where Q is the source strength and a. is the standard deviation of 

i. 1 yv ‘ — _ —-4 ^ 4— 4*" -1 M 4- ■« 4 4 4" -I 1 “4 ^ ^ J O 4" “4 « 4 T T /A ^ O 4“ A 4 T 


the concentration distribution in the 
dosage term is given by 


direction; whereas, the peak 


2tt u a a 

y z 

where u is the mean wind speed. The second difference between these 
models is that the concentration contains a modular alongwind term to 
account for downstream temporal effects not considered in the dosage 
model. The alongwind term affords an exponential decay in concentration 
as a function of: cloud transit time, concentration distribution and 
the mean wind speed. 


The lateral term - which is common to both models - is 
another exponential decay term, and is a function of the Gaussian 
spreading rate and the distance laterally from the mean wind azimuth. 
The vertical term - again common to both models - is a rather complex 
decay function since it contains a multiple reflection term for the 
point source which stops the vertical cloud development at the top of 
mixing layer and eventually changes the form of the vertical concen- 
tration distribution from Gaussian to rectangular. The last modular 
in both models is the depletion term. This term accounts for the loss 
of material by simple decay processes, precipitation scavenging or 
gravitational settling. 
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Each of the modular terms In these two relations is 
heavily dependent upon empirical parameters. For this reason, it is 
necessary to perform reliable investigations of as many rocket launches 
as possible and to perform controlled investigations of rocket exhausts 
experiments to insure the maximum reliability in these empirical 
parameters. This will in turn afford a more accurate dispersion 
description from these models. 

c. The Description of the Models in the NASA/MSFC Multi- 
layer Diffusion Model 

The normal launch environment will usually involve an 
atmospheric structure comprised of several horizontal meteorological 
layers with distinctive wind velocity, temperature, and humidity regimes 
between the surface and a 5 kilometer altitude. Large horizontal 
spatial variation in these meteorological parameters may also occur in 
the surface layer as a consequence of changes in terrain or land-water 
interfaces, which is accounted for by the diffusion model. The general 
diffusion model for the concentration (Equation 1) and the dosage 
(Equation 2) assumes an expanding volume about a moving point of 
reference in a homogeneous environment. 

To overcome the obvious shortcomings of the general 
diffusion model but to stay within the established bounds of classical 
fluid mechanics [2], a multiple layer concept is introduced to cope 
with the vertical and horizontal atmospheric gradients. Here, the 
general diffusion model is applied to individual horizontal layers in 
which the meteorological structure is reasonably homogenous and indepen- 
dent of the neighboring layers. These layers have boundaries which are 
placed at points of major discontinuities in the vertical profiles of 
wind velocity, temperature, and humidity. Since the Multilayer Diffusion 
Model has imposed the general restriction of layer independence (no 
flux of particles or gases entering or leaving an individual layer), 
special provision must be made for spatial changes in the horizontal 
meteorology and for gravitational settling or precipitation scavenging. 

In addition, the type of source within a layer must be considered; 
that is, whether there is a ground cloud source or a plume cloud 
source (see definitions). 

The NASA/MSFC Multilayer Diffusion Model has six models 
(Figure 1) which account for three categories of dispersive constraints: 
The source distribution, the environmental effects and the depositional 
effects. This flexibility is required to deal with the stages of the 
development of the exhaust cloud and the complex potentially varying 
meteorological conditions. These models can be used alone to 
describe all the environmental layers or in superimposed combinations 
where variations in layer meteorology require different modeling. For 
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METEOROLOGICAL INPUTS 


SOURCE INPUTS 


CONCENTRATION, DOSAGE, AND DEPOSITION MODELS 

I. SOURCE DISTRIBUTION 

1. ELLIPTIC-CYLINDRICAL SOURCE EXTENDS VERTICALLY 
THROUGH ENTIRE DEPTH OF LAYER AND TURBULENT MIXING 
IS OCCURRING 

3. ELLIPSOIDAL SOURCE DOES NOT EXTEND VERTICALLY 
THROUGH ENTIRE DEPTH OF LAYER 

II. ENVIRONMENTAL EFFECTS 

2. SOURCE EXTENDS VERTICALLY THROUGH ENTIRE DEPTH OF 
LAYER AND TURBULENT MIXING IS NOT OCCURRING 

4. FULL TRANSITION MODEL FOR STEP-CHANGE IN LAYER 
STRUCTURE 

III. DEPOSITION EFFECTS 

5. PRECIPITATION SCAVENGING 

6 . GRAVITATIONAL SETTLING 


LOGIC SECTION 1 
CALCULATES DOSAGE AND CONCEN- 
TRATION PATTERNS AND SURFACE 
DEPOSITION DUE TO PRECIPITATION 
SCAVENGING: MODELS 1, 2, 3, 4, 5 


LOGIC SECTION 5 
CALCULATES DEPOSITION DUE 
TO GRAVITATIONAL SETTLING: 
MODEL 6 


LOGIC SECTION 2 
CALCULATES PEAK 
DOSAGE AND PEAK 
CONCENTRATION: 
MODELS 1, 2, 3 



LOGIC SECTION 4 
CALCULATES ISOPLETHS 
OF DOSAGE AND CONCEN- 
TRATION IN THE y-z 
PLANE: MODELS 1, 2, 3 


- 

r 



LOGIC SECTION 3 1 

CALCULATES ISOPLETHS OF DOSAGE 

AND CONCENTRATION IN THE x-y 

PLANE: MODELS 1, 2, 3 

1 ' 


FIGURE 1. 


j jOUTPUT LISTING J 

BLOCK DIAGRAM OF THE COMPUTER PROGRAM FOR THE NASA/MSFC 
MULTILAYER DIFFUSION MODEL. 


13 












the introductory overview, however, these combinations will not be 
considered. The primary output of all models is a mapping of the 
regimes of the concentration and dosage isopleths and centerline profiles 
for concentration and dosage. 

The fundamental category of dispersive constraints is the 
source distribution . The two distributions are: 

1. The e lliptic-cyclindrical source which assumes a 
two dimensional Gaussian distribution in the x-y plane and a uniform 
distribution in the vertical direction. 

2. The ellipsoidal source which assumes a three- 
dimensional Gaussian distribution. 

Model 1 is for the elliptic-cylindrical source whose 
vertical expansion is constrainted by the layer boundaries -- thus has 
only a two dimensional expansion in the horizontal plane due to tur- 
bulence mixing. This model is normally used to describe the rocket* s 
inflight plume cloud. 

Model 3 is for the ellipsoidal source and is assumed to 
expand in all three dimensions as the effluents are propagated down- 
stream. When the e llipsoidica 1 source reaches the top of the mixing 
layer, the distribution of the constituents is reflected back into 
the expanding vertical distribution. On the other hand, that fraction 
not lost in surface deposition is also reflected back in a similar 
manner. After sufficient mixing, the ellipsoidal distrituion becomes 
a elliptic-cylindrical distribution (Model 1). While Model 3 is 
normally used to describe the dispersion of the rocket * s ground cloud, 
it could be used to model upper air explosions. The formulation for 
Model 3 has been provided in Appendix II. 

The second category is environmental effects . The two 

effects are: 


1. No turbulence mixing in the upper atmosphere. 

2. Changes in meteorological conditions as the con- 
stituents are transported downstream. 

Model 2 is the same as Model 1 except it is assumed that 
there is no turbulent mixing. This implies that the exhaust material 
just meanders along the layer without dispersing. While the Model 2 
is not generally used, movies of rocket firings clearly show that under 
some special meteorological conditions this model is required. While 
the Multilayer Diffusion Model is general in applicability, it is 
specific in meteorological parameters and launch description. 
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Model 4 updates the diffusion model with changes in 
meteorological conditions and structure which can occur as the con- 
stituents propagate downstream. This model assumes that the vertical 
concentration of material has become uniform throughout each layer when 
a step-change in the meteorological conditions is introduced, resulting 
in the destruction of the original layer boundaries and the formation 
of new layer boundaries. The concentration fields which exist at this 
time are treated as new sources. In those new layers which now comprise 
more than one old layer, the old concentration is mapped as two indepen- 
dent concentration sources and then superimposed for the resulting 
concentration and dosage mappings. 

The third category of dispersive constraints includes the 
deposition due to: 

1. Precipitation scavenging 

2. Gravitational settling 

Model 5 accounts for precipitation scavenging. An example 
of where Model 5 must be used is in solid rocket launches during the 
occurrence of rain, because the HC1 will be scavenged by the rain. 

Model 6 describes the ground deposition due to gravitational settling 
of particles or droplets. Wind shears are incorporated in this model 
to account for the effect of the settling velocity of the particulate 
matter. There are two forms for the source in this model; namely: 

1. The source that extends vertically through the 
entire layer with a uniform distribution — this is the same source 
model as used with Models 1 and 2, and 

2. A volume source in the K-th layer -- this is the 
same source model as used with Model 3. 

Model 6 is very important in the analysis of the settling 
of Al^Og particles released in solid rocket firings. 

The treatment of cold spills and fuel leaks that occur 
near ground level requires a continuous source, but the models that 
have been considered so far are for discrete sources; therefore, the 
models must be adapted for the use in predicting concentration-dosage 
levels downwind from continuous sources. 

The layer of the environment influenced by the ground-level 
spills and leaks can be treated as homogeneous; therefore, the general 
formula for concentration and dosage (Equations 1 and 2) presented in 
the initial discussion would be applicable if spills and leaks are 
treated as continuous sources. To visualize this adaption for these 
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formulas, assume a source cloud with a concentration distribution 
that implies a given dosage at a point for this cloud; that is, the 
dosage per event. If there are a number of similar clouds, discretely 
spaced, then for each cloud we obtain a dosage for each cloud whose 
sum corresponds to the total dosage for the entire event. 

In the limit as the spacing between clouds approaches zero 
and the number of clouds becomes large, the discrete source approaches 
a continuous source whose concentration is the point dosage per unit 
time. The relation for the continuous concentration, which follows 
directly from this argument, is 

Concentration = (Peak Concentration) x (Lateral Term) 

x (Vertical Term x (Depletion Term) 

X - > ( 3 ) 

^peak 2tt u a a L 

^ y z R 

where Q is the source strength, u is the mean wind speed, a. is the 
concentration distribution in the i-th direction and t is the 

"D 

release time. 


The Lateral Term, Vertical Term, and the subset of 
equations defining a and cr are the same as for the point source 

y 2 

dosage (Equation 2). The continuous dosage is then the continuous 
concentration times the release time. 

In summary, the Multilayer Diffusion Model is composed of 
six submodels. Models 1 and 3 are designed to distinguish between the 
two sources of toxic cloud formation -- the ground cloud during the 
initial launch phase (Model 3) and the plume cloud after the initial 
launch phase (Model 1). From the stand point of environmental impact, 
the description of the fields of the ground deposition of materials 
from the ground cloud is of primary significance -- this description 
is afforded by Model 3. Generally, this model is employed in the 
surface layer, but can be employed in any layer where the source does 
not extend through the entire layer. 

Model 2 was designed to account for a lack of turbulent 
mixing which can occur in the upper atmosphere. Model 4 is employed 
when a change in meteorological condition occurs during the downstream 
transport of the cloud. In the event of rain, the precipitation 
scavenging — both of gases and particles -- can be accounted fc^f in 
Model 5. The fallout of particulate matter on the ground is the domain 
of Model 6. These six submodels form the basic set of equations which 
are available to treat the diffusion problem. To model a specific 
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launch of a vehicle, it is necessary to blend these equations together 
and adjust the model parameters to the specific meteorological conditions 
of the launch, to the specific terrain around the launch site and to 
the specific vehicle being launched; thus the degree of complexity 
in the diffusion model. 
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SECTION III. PREDICTIONS OF THE DISPERSIVE TRANSPORT OF SATURN V 
ENGINE EXHAUST EFFLUENTS IN THE TROPOSPHERE (APRIL 16, 1972) 


Concentration mappings of the exhaust effluents from 
the launch of the Apollo 16 mission (Apollo Saturn 511) at 1754 Univer- 
sal Time (UT) (1254 EST) from Kennedy Space Center obtained from the 
NASA/MSFC Multilayer Diffusion Model are presented in this section. 

Since the calculation of the exhaust effluent dispersions presented in 
this report are limited to altitudes below 18 km, only the exhaust 
effluents from the first stage (S-1C) of this Saturn V have an effect 
upon the concentration mapping and the effluents from the remaining 
stages can be neglected here. 

The two important input parameters in the Multilayer Diffusion 
Model are the vehicle fuel properties and the meteorological conditions; 
these input parameters will be disaussed as an introduction to the 
actual model parameters employed and the predictions obtained. 

a. Saturn V (S-1C) Engine Exhaust Effluents 

It has been established theoretically and from observa- 
tions that low altitude clouds of engine exhaust and entrained gases 
are formed during the launch of the S-1C booster [3], In the case of the 
S-1C vehicle, these clouds regularly rise and stabilize at an altitude 
of about 2 km, depending upon meteorological conditions. The engine 
exhaust effluents contained in this cloud are produced during the period 
from rocket ignition to the time that it ascends to the altitude as 
which the exhaust ground cloud stabilizes. For a cloud stabilization 
altitude of two km, the ground cloud contains the exhaust effluents 
produced during approximately the first 43 seconds of engine ignition, 
or about 6.524 x 1CT kg (Table I). The composition of these liquid 
engine exhaust effluents at engine exit and after reaction with the 
atmosphere consists primarily (99.9 percent, see Table II) of water 
vapor (H 2 O) , carbon monoxide (CO) and carbon dioxide (CO 2 ) • Of these 
gases, the only two that are considered toxic are carbon monoxide and carbon 
dioxide (see Appendix II, Toxicity Criteria). The stabilized ground cloud 
from Apollo 16 contained about 2.85 x 10^ kg of carbon monoxide [4]. 

g The single-start F-l engine, with a fixed thrust of 6.8 x 

10 N, uses liquid oxygen as the oxidizer and RP-1 as the fuel. Water, 
hydrocarbons, and other carbon compounds make up the majo^ bulk of the 
engine exhaust. Eight minutes before lift-off, 5.02 x 10 kg/sec 
(8,000 gallons/minute) of water are sprayed into the trench (Figure 2). 

At 2 minutes before lift-off, 1.26 x 10^ kg/sec (20,000 gallons /minute) 
of water are sprayed onto the mobile launcher, launch pad, and into the 
trench for cooling purposes. Six seconds after ignition, the F-l 
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TABLE I 


EXHAUST MATERIALS EMITTED AS A FUNCTION OF ALTITUDE FROM THE S-1C STAGE 


— ______ — - 

Altitude Range 

Range Time at Top of Layer 

Exhaust Material Emitted 

(kin) 

(sec) 

(kilograms) 

0-2 

36** 

5.460 x 10 5 

2-4 

49 

1.690 x 10 5 

4-6 

58 

1.170 x 10 5 

6-8 

66 

1.040 x 10 5 

8-10 

72 

7.800 x 10 4 

10-12 

78 

7.800 x 10 4 

L2-14 

83 

6.500 x 10 4 

14-16 

88 

6.400 x 10 4 

16-18 

92 

5.200 x 10 4 

18-20 

96 

5.200 x 10 4 

20-22 

100 

5.200 x 10 4 

22-24 

104 

5.200 x 10 4 

24-26 

107 

3.900 x 10 4 

26-28 

110 

3.900 x 10 4 

28-30 

113 

3.900 x 10 4 

30-40 

128 

1.950 x 10 5 

40-50* 

140 

1.456 x 10 5 

50-60 

151 

1.144 x 10 5 

60-70 

160 

9.360 x 10 4 


*Layer in which IBECO occurs; only four engines burn after T + 135.96 
seconds. 


**Total effluent calculations based on 43 seconds in that ignition 
occurs at T - 7 seconds. 
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TABLE II 


F-l EXHAUST COMPOSITION 



Component 

Weight Percent 
at Exit 

Estimated Weight 
Percent After Reaction 
w/Atmosphere * 

H 2 

1.13 


OH 

.22 


H 2° 

26.6 

27.5 

0 

.005 


°2 

.011 


CO 

42.2 

43.6 

C°2 

27.8 

28.8 

CHO 

.002 


H 

.022 


HC 

1.60 


Particulates 

.0575 

.0575 

NO ! 

x 

0 

.0005 


Hydrocarbons and particulates result primarily from turbine exhaust. 
Hydrocarbons are expected to burn when mixed with primary flow and 
atmosphere. 

Primary exhaust products are based on theoretical data. Imperfect mix- 
ing is the primary mechanism which might cause the measured composition 
to differ from theoretical. 

Measured turbine exhaust composition agrees very well with predicted 
composition but no data have been located to confirm the expected 
reaction with primary flow and atmospheric oxygen. 

*These estimates were made by S&E-AERO-AT , Marshall Space Flight 
Center, Alabama, National Aeronautics and Space Administration. 
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FIGURE 2. CUTAWAY VIEW OF LAUNCH COMPLEX 39A 



engines produce about 1.34 x 10 calories of heat per second. However, 
after evaporative cooling from the water spray during initial launch, 
only about 1.27 x 10^ calories of heat per second are emitted [4]. 

The exhaust cloud from the Apollo 16 launch was monitored 
photographically at two-second intervals to provide empirical data to 
evaluate cloud rise calculations from the Multilayer Diffusion Model. 

The cloud rise and growth rate data were extracted from 
the sequential set of photographs by using scaling factors to account 
for the distance between the camera and the clouds, for the diameter of 
the top of the clouds and for time (Figures 3 and 4), 

The predicted maximum cloud rise and the fitted curves of 
the observed rates of the bimodal cloud for Apollo 16 launch are shown 
in Figure 3. The greater rise rate of the right plume is due to 
deflection by the 21-degree inclination of the south end of the concrete 
flame trench. However, the two clouds eventually rise to an equal 
height and unite at the time they become stable with the ambient air. 

b. Atmospheric Conditions at the Time of Launch of Apollo 16 

Predictions of the concentration and movement of rocket 
motor effluents require inputs of several meteorological parameters 
varying in complexity from relatively straight forward ones such as 
potential temperature and wind velocity to rather complex ones such as 
standard deviations of the wind azimuth and elevation angles at all 
altitudes of interest. Selection of meteorological inputs begins with 
the assignment of layer boundaries based upon the vertical profiles of 
wind, temperature, and humidity within this volume. 

The sources for meteorological information at Kennedy 
Space Center at the time of launch of Apollo 16 included the *Jimsphere 
and radiosonde measurement systems, United States Air Force weather 
reconnainssance aircraft operating in the area, and the NASA 150-meter 
Ground Wind Tower on Merritt Island. Other sources of atmospheric 
data included synoptic weather maps and upper altitude atmospheric 
pressure and wind charts. Meteorologists from MSFC started analyses 
of meteorological data about two days before the launch of Apollo 16 to 
make preliminary plans for aircraft sampling profiles. These analyses 
continued up to the time of launch. 

The atmospheric conditions in the vicinity of Kennedy 
Space Center Launch Complex 39 (LC39A) at the time of launch on April 16 
were influenced primarily by dry stable air associated with an anti- 
cyclonic system centered off the east coast of Florida. Locally, there 
were only a few scattered cumulus clouds with bases at about two km. A 
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FIGURE 3. RISE RATE OF BIMODAL EXHAUST GROUND CLOUD FOR THE LAUNCH OF APOLLO 16 ON APRIL 16, 
1972 (1754UT) AT KENNEDY SPACE CENTER, FLORIDA 


PREDICTED BY CLOUD 



FIGURE 4. HORIZONTAL GROWTH OF GROUND CLOUD FOR THE LAUNCH OF APOLLO 16 
ON APRIL 16, 1972 (1754UT) AT KENNEDY SPACE CENTER, FLORIDA 



sea breeze had penetrated inland from the southeast across much of Cape 
Kennedy but had not reached as far as LC 39A. The wind was southwest at 
almost all altitudes from the surface to an altitude of 15 km (Figure 5). 
These persistent southwest winds were responsible for delaying and 
reducing the effect of the sea breeze. Wind speeds ranged from 5.0 m/sec 
at the surface to 26 . 1 m/sec at 11.9 km altitude. In addition to wind 
data, a crucial parameter influencing the altitude at which the ground 
exhaust cloud will stabilize is atmospheric stability. At the time of 
the Apollo 16 launch, a subsidence inversion existed over KSC , with 
a base of about 2000 m. This inversion had a profound effect on the 
ground exhaust cloud dynamics, and subsequently on the concentration of 
rocket engine effluents both at the ground cloud and ground surface 
levels . 


The southwesterly wind pattern developed over Florida as 
a consequence of a low pressure trough moving eastward from the 
Mississippi Valley. As this trough shifted to the higher pressure 
existing over Florida, the pressure gradient increased through the 
troposphere resulting in fairly strong southwesterly winds. A rapidly 
moving cold front associated with the trough aloft was located over 
extreme northwest Florida at launch time, but the rain shower activity 
produced by this weather system was confined to the immediate vicinity 
of the front; therefore, the front had no direct effect on the local 
weather conditions. However, the influence of the weather system 
dominated the wind profile, and this factor was critical in making 
predictions of the ground exhaust cloud movement. 

c. Meteorological and Source Inputs to the Model 

Based on the vehicle used, in this case the Saturn V, and 
the meteorological conditions at the launch time, a set of values for 
meteorological and source model input parameters for the MSFC/NASA 
Multilayer Diffusion Model are defined. Time-lapse photographs are used 
primarily as an empirical standard comparison for the observed and predicted 
exhaust cloud rise and volume. A detailed list of the general input 
parameters required by the model is given in Appendix III. The specific 
values that were employed in the model to obtain the concentration- 
dosage mapping for the launch of the Saturn V vehicle in the Apollo 16 
mission are given. 

Calculations of the maximum height (z ) of rise of the 
Apollo 16 engine exhaust ground cloud, to define the location and size 
of the source, were made using rawinsonde soundings at 6 hours 39 
minutes before launch (1115 UT) to forecast and at 10 minutes after 
launch (1804 UT) to predict the environmental dispersion of these 
exhaust effluents. The values of the ambient temperature (T) , the^ 
potential temperature gradient (30/Sz)> and the mean wind speed (u) 
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for these two soundings are given in Table III. This positive 
potential temperature gradient implies that a stable cloud rise relation 
is required (Equation A-2 or A-4) . The relatively slow rise of the 
Saturn V ground Cloud dictates the use of the continuous cloud rise 
relation for the stable source. The initial vertical plume velocity 
(w ), the initial radius (r ), the entrainment coefficient (Y ) (Table III) 
are required in the cloud rise calculation and are based on the charact- 
ertistics of the Saturn V vehicle [4 ] . The maximum cloud rise calculated 
(Equation A-4) for the 1804 UT sounding is 2000 meters. This analytical 
prediction of the height in Figure 3 is compared to the observed cloud 
rise heights (solid curve) obtained using time-lapse photographs [5 ]• 
(While Figure 3 only gives a single analytically predicted cloud rise 
height, a predicted cloud rise curve can be analytically generated 
with Equation A-4.) 


The height of cloud stabilization is then employed in the 
empirical prediction of the cloud diameter (Figure 4) based on the MSFG library 
of cloud growth observations. For comparison of the predicted with 
observed cloud diameter, the photographically observed cloud growth 
as a function of time was made. 


The ground cloud became stable with the ambient air 150 
seconds after launch and had a maximum height of about 2000 meters. At 
this time, the cloud contained 5.460x10^ kilograms of the rocket engine 
effluents from approximately the first 43 seconds of engine operation. 
Calculations of the parts per million (ppm) of carbon monoxide in the 
stabilized cloud indicated about 130 ppm if uniform mixing occurred; 
however, in the application of the MSFC/NASA Multilayer Diffusion Model 
it is assumed the exhaust effluents in the cloud are Gaussian distributed 
in the plane of the horizon, resulting in a higher centerline concen- 
tration of carbon monoxide than the average value just given. Photo- 
graphs of the cloud at 6.0 minutes after launch revealed that the ground 
cloud had assumed a shape modeled by oblate spheriodal (Figure 6). The major 
axes was 2000 meters in the x and y directions and the minor axis was 
800 meters in the z direction (the x-axis is in the alongwind direction; 
y is in the crosswind direction; and z is in the vertical direction). 

A detailed examination of the appropriate air temperature 
and dew point (Figure 7), and wind speed and direction (Figure 8) at 
the time of the launch of Apollo 16 indicated the atmosphere should be 
initially assumed to consist of nine layers^ The top of seven of these 
layers are defined by the meteorological parameters indicated in these 
figures, and the other two layers are defined by cloud dimensions. Initially, 
the MSFC Multilayer Atmospheric Diffusion Model 1 (source extends vertically 
through the entire depth of the layer and vertical mixing is occurring) was 
applied to initiate the distributions. After stabilization of the ground 
exhaust cloud, it was assumed there was layer breakdown between the lowest 
six layers, and therefore turbulent mixing would occur in the first 2000 
meters above the launch pad. Model 4 (full transition model for step-change in 
layer structure), with appropriate layer breakdown formulas was used to simulate 
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FIGURE 6 


IOOO 600 200 0 200 600 

CROSSWIND CLOUD DIMENSION (meters) 


. CROSSWIND AND VERTICAL CLOUD DIMENSIONS ASSIGNED TO 
THE GROUND CLOUD FORMED FROM THE LAUNCH OF APOLLO 16 










• - POINTS WHICH DEFINE LAYER BOUNDARIES 



FIGURE 7. AIR TEMPERATURE AND DEW-POINT TEMPERATURE PROFILES FROM 
THE 1804 UT RAWINSONDE ON APRIL 16, 1972 AT CAPE KENNEDY 


30 



• - POINTS WHICH DEFINE LAYER BOUNDARIES 



WIND DIRECTION (degrees) 


FIGURE 8. WIND SPEED AND DIRECTION PROFILES FROM THE 1804 UT 
RAWINSONDE ON APRIL 16, 1972 AT CAPE KENNEDY (SOLID 
LINES). DASHED LINES REPRESENT ADJUSTED PROFILES 
BASED ON INFORMATION FROM THE 150 -METER GROUND WIND 
TOWER 
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the effects of this turbulent mixing of the six lowest layers. Between 
altitudes of 2000 meters and 4000 meters (altitude of the next inversion 
layer), Model 3 (vertical extend of the source is less than the depth 
of the layer) was assumed. This resulted in trapping of all engine 
effluents of the ground exhaust cloud that were above an altitude of 
2000 meters within the altitude region 2000 to 4000 m. The portion of 
the ground exhaust cloud contained between altitudes of 1800 and 2000 
meters was assumed to mix and diffuse downward to the ground. 

Source and meteorological input parameters for the MSFC/ 

NASA Multilayer Diffusion Model are defined in Appendix III (Tables 
A- 1 and A-III) . The values of the meteorological input in Table IV 
are for all nine layers prior to the breakdown of the layers. Values 
for the source dimensions and source strengths in a respective layer 
were obtained relative to the stabilized ground cloud (Figure 6) 
utilizing the relations in Appendix HI (Equation A-39, A-40, and A-41) . 
The resulting source parameter for the launch of Apollo 16 are given 
in Table V. These meteorological and source input parameters (Tables 

IV and V) are the basis of the effluent dispersion results for the 
launch of Apollo 16 that were obtained with the NASA/MSFC Multilayer 
Diffusion Model. 

d. Predictions of the Concentration Contours from Exhaust 

Effluents for Apollo 16 by the NASA/MSFC Multilayer 

Diffusion Model 

The following predictions for the dispersion of the Saturn 

V exhaust effluents from the Apollo 16 mission will focus on the 
dispersion of carbon monoxide and carbon dioxide because these are the 
only two constituents that are potentially toxic and have a potential 
impact on the environment air quality in the troposphere. Carbon 
monoxide is selected as the trace gas in the isopleths and centerline 
concentration profiles that are discussed because the concentrations of 
carbon monoxide are 2.4 times stronger than those of carbon dioxide 
and because the toxicity levels for carbon monoxide are much less than 
those for carbon dioxide (Appendix IV); that is, if there is not an 
environmental hazard from carbon monoxide during the launch of a Saturn V 
vehicle, then it can be assumed that the other constituents in the exhaust 
effluents will not afford an environmental hazard. 

To ascertain the surface concentrations of carbon monoxide 
a nine layer atmosphere structure (Figures 7 and 8) was employed, based 
on the atmospheric thermodynamic and kinematic soundings at launch time, 
initially in Model 1 of the NASA/MSFC Multilayer Diffusion Model (as 
discussed in the last section) and run for one second. Then Model 4 
was invoked to merge the first six layers. The top of this new surface 
layer was dictated by the temperature inversion at 2 kilometers. The 
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TABLE IV 


METEOROLOGICAL MODEL INPUTS FOR THE NINE LAYERS 







Layer 





Parameter 

Units 












m 

2 

i. . . 


n 

5 

6 

7 

8 

9 

u 

R 

m sec ^ 

5.24 









a 

BK 

-i 

in sec 

5.0 

9.3 

9.6 

9.4 

10.1 

9.9 

11.1 

10.5 

8.1 

u 

TK 

-1 

m sec 

9.3 

9.6 

9.4 

10.1 

9.9 

11.1 

10.5 

8.1 

6.7 

e 

BK 

o 

K 

301.0 

301.9 

301.7 

301.8 

301.9 

301.8 

302.1 

305.8 

308.3 

e 

TK 

o 

K 

301.9 

301.7 

301.8 

301.9 

301.8 

302.1 

305.8 

308.3 

309.5 

7 {T } 
AR oK 

deg 

8.66 









7 {t } 

ABK OK 

deg 

9.08 

4.92 

4.54 

4.34 

5.19 

4.14 

4.09 

2.0 

1.0 

a {t } 
ATK L or 

deg 

4.92 

4.54 

4.34 

4.19 

4.14 

4.09 

2.0 

1.0 

1.0 

T 

oK 

sec 

600 

600 

600 

600 

600 

600 

600 

600 

600 

y 

K 


1 

■! 

1 

1 

1 

1 

1 

1 

1 

1 


33 













TABLE V. 


SOURCE INPUTS FOR STABILIZED CLOUD 
















resulting isopleths are shown in Figure 9. These isopleths represent 
the spatial mapping at the bottom of the original six layers. This 
spatial selection represents the significant locations (defined by 
meteorological conditions) required in an environmental inventory of 
the impact of the rocket exhaust effluents. Of obvious importance are 
the surface (2 meter level) isopleths, which reveal that the maximum 
concentration of carbon monoxide was less than 2 ppm from the launch of 
Apollo 16. This is well below the public exposure level of 25 ppm in 
an one-hour period recommended by the Environmental Protection Agency 
(Appendix IV, Table A-III) . The maximum centerline downwind concentrations 
of CO and CO at the bottom of the layers between the source center 
(2.2 kilometers) and the surface are shown in Table'VI. Turbulent 
mixing and diffusion processes diluted the constituents rapidly, as can 
be observed from the small quantities that reached ground level. To 
convert concentrations of CO (as indicated by the isopleths in Figure 9) 
to CO 2 , it is necessary to multiply by a conversion factor of .4203. 

The factor was determined from the following relationship: 


28.8% CO 

ppm (CO) x x molecular WT of CO = pp m (CO ) (4) 

43.6% CO molecular WT of 2 

The 28.8% CO 2 and 43.6% CO values in the above equation represent the 

estimated percentage by weight of these Saturn engine exhaust constituents 
after reaction with the atmosphere. 

Four downwind maximum centerline concentrations of carbon 
monoxide were calculated, ranging in alitutde from 2.2 kilometers, 
which was the center of the stabilized ground cloud, to 17.5 kilometers. 
The upper three altitudes (3.5 km, 8.5 km, and 17.5 km) were selected 
to correspond to the alitutdes at which Air Force aircraft obtained 
samples. The ground cloud (2.2 km) was treated as an oblate spheroid 
(this is the original ground cloud) utilizing Model 3. In the remaining 
three levels, the source is the plume cloud of uniform cylindrical con- 
centration. The cylinders are assumed to have a height of 1.0 kilometer 
and a diameter of 400 meters after a cloud stabilization time of 30 
seconds (Figure 10). The height of these layers affect the concentration 
because Model 1, which assumes a uniform vertical concentration, is used. 
The concentrations downwind are given in Figure 11 for these four 
altitudes. This figure shows that the downwind concentrations at 8.5 
kilometers and 17.5 kilometers altitude are greater than those at 3.5 
kilometers beyond a downwind distance of 20 kilometers; showing clearly 
the support transport effect of the isothermal layer at an altitude of 
about 4,0 kilometers. The reason that the concentrations decreases as 
a function of increasing altitude is that the vehicle gains speed as its 
altitude increases and therefore spedds less time in each layer as it 
ascends . 
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GROUND CLOUD DATA 



*Height at center of stabilized cloud 

^Concentrations are high at center of cloud because of gauss ian assumptions 


















Altitude (km) 
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In summary, these calculations indicated maximum ground 
level concentrations of CO of less than 2 ppm for the Apollo 16 launch, 
well below the 25 ppm maximum allowable one- hour concentrations allowed 
by the Environmental Protection Agency, The maximum calculated con- 
centration of CO in the ground exhaust cloud at an altitude of 2200 
meters was 295 ppm, rapidly diminishing to 23 ppm by turbulent mixing 
and diffusion processes at a downwind distance of about 30 km. 


SECTION IV. PREDICTIONS OF THE DISPERSIVE TRANSPORT OF SCOUT -ALGOL III 
MOTOR EXHAUST EFFLUENTS IN THE TROPOSPHERE (AUGUST 13, 1972) 


Concentration mappings of the exhaust effluents from the launch 
of the Scout-Algol III vehicle at 1510 Universal Time (UT) (1110 EDT) 
from Wallops Island, Virginia, obtained from the NAS A/ MS F C Multilayer 
Diffusion Model are presented in this section. The calculation of the 
exhaust effluent dispersion presented in this report are limited to 
altitudes below 500 meters because the exhaust effluents released by 
the Scout vehicle are relatively small; therefore, only the ground cloud 
affords a sufficiently intensity source of effluents for a meaningful 
analysis a t the surface. In contrast to the Saturn V exhaust effluents, 
the Scout exhaust effluents contain more potentially toxic constituents and, 
obviously, because of the size of the Scout, the concentrations of 
these efflunets are much less. 


As it was illustrated in the last section, the two important 
parametric inputs in the MSFC/NASA Multilayer Diffusion Model are the 
vehicle fuel properties and the meteorological conditions; these Scout- 
Algol input parameters will be discussed as an introduction to the actual 
model parameters employed and the predictions obtained. 


a. Algol III Motor Exhaust Effluents 


It has been established theoretically and from observa- 
tions that low altitude clouds of motor exhaust and entrained gases are 
formed during the launch of the Algol III booster. In the case of the Scout- 
Algol III vehicle, these clouds regularly rise and stabilize at an altitude 
of about 400 meters (based on empirical observations) depending upon 
meteorological conditions. The motor exhaust effluents contained in 
this cloud are produced during the period from rocket ignition to the 
time that it ascends to the altitude at which the exhaust ground cloud 
stabilizes. For a cloud stabilization altitude of 400 meters, the 
ground cloud contains the exhaust effluents produced during the first 
eight seconds of motor ignition, or about 1592 kg of material 
(Table VII). The composition of these solid rocket motor exhaust 
effluents consists primarily (98.3 percent, see Table VIII) of, in order of 
descending concentration, alumina (AI 2 O 3 ), carbon monoxide (CO), hydrogen 
chloride (HC1) , nitrogen (^) > water vapor (1^0) , hydrogen (H 2 ), and 
carbon dioxide- (CO 2 ) . Four are potentially toxic (see Appendix IV). The 
stabilized ground cloud from the Scout-Algol III contained about 511 kg 
of alumina, 436 kg of carbon monoxide, 328 kg of hydrogen chloride, and 
only 33 kg of carbon dioxide. 


The rocket plume from the Scout-Algol III before lift-off 
impinges on a steel plate surrounded by concrete in contrast to the flame 
trench with water that is employed for Saturn launches. The resulting 
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TABLE VIII. 


FUEL PROPERTIES OF THE SCOUT-ALGOL III FIRST-STAGE MOTOR 


Exhaust Products at 

Exit (percent by weight): 


Alundum 

- ai 2 o 3 * ** 

32.1 

Carbon Monoxide 

- CO** 

27.4 

Hydrogen Chloride 

- HC1** 

20.6 

Nitrogen 

“ N 2 

7.2 

Water Vapor 

-HO 

2 

6.4 

Hydrogen 

- H 

2 

2.5 

Carbon Dioxide 

- co 2 ** 

2.1 


Other 

1.7 


-1 5 

Fuel Expenditure Rate (g sec ) 1.9414 x 10 * 

Fuel Heat Content (cal g" 1 ) 979.4 

Total Burn Time (sec) 78 


*Exact values are given in Table VII - this value is the average over 
the initial 13.3 seconds of ascent. 

^Potentially toxic (see Appendix IV) 
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Scout-Algol III plume is cylindrically symmetrical rather than a bimodal 
Saturn V plume. Another difference between the Scout and Saturn vehicle 
launches is that the Scout lift-off is milliseconds after ignition, 
rather than the 7 seconds associated with the Saturn. As a result of 
these factors, the vast amount of entrained water found in the Saturn V 
ground cloud is not present in the Scout ground cloud. 

The exhaust cloud from the Scout-Algol III launch was 
monitored photographically at five-second intervals to provide empirical 
data to evaluate cloud rise calculations from the Multilayer Diffusion 
Model. Two cameras were employed in an orthogonal configuration west 
and south of the launch area, pointing toward the launch area. 

The cloud rise and growth rate data were extracted 
from the sequential set of photographs by using scaling factors to 
account for the distance between the camera and the clouds, for the 
diameter of the top of the clouds and for time. 

b. Atmospheric Conditions at the Time of Launch of 

the Scout-Algol III 

Predictions of the concentration and movement of rocket 
motor effluents require inputs of several meteorological parameters 
varying in complexity from relatively straight forward ones such as 
potential temperature and wind velocity to rather complex ones such as 
standard deviations of the wind azimuth and elevation angles at all 
altitudes of interest. Selection of meteorological inputs begins with 
the assignment of layer boundaries based upon the vertical profiles of 
wind, temperature, and humidity within this volume. 

The sources for meteorological information at Wallops 
Island at the time launch of the Scout included rawinsonde measurements, 
and the two NASA 76-meter Ground Wind Towers. Other sources of atmospheric 
data included synoptic weather maps and upper altitude atmospheric 
pressure and wind charts. Meteorologists from MSFC started analyses of 
meteorological data about two days before the launch of the Scout-Algol III 
vehicle to make preliminary plans for aircraft sampling profiles. These 
analyses continued up to the time of launch. 

The atmospheric conditions in the vicinity of the 
Wallops Island launch site at 1510 UT on August 13, 1972 were dominated 
by partially cloudy skies, light winds and a surface temperature of 25°C. 

No strong surface pressure systems were in the immediate area, although 
a weak cold front over northern Virginia was moving toward the launch 
site. Precipitation was not evident in the Wallops area at launch time 
and cloud layers were primarily high (4km to 6km). This frontal system 
was in conjunction with a short wave of low pressure aloft moving through 
the long wave pattern established over the east coast of the United States 
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and the western Atlantic Ocean. The surface wind was southerly at 3,9 m/sec. 
The rawinsonde sounding (Figures 12 & 13) taken 40 minutes prior to 
launch (1430 UT) , indicated a baroclinic spiraling of the wind direction 
such that at 2 kilometers the wind direction had shifted to a westerly 
direction. The wind speed, however, decreased with altitude to a near 
calm at the base of the elevated inversion (1.2 kilometers): at the top 
of this inversion (1.5 kilometers) it increased with altitude to about 
3 m/sec at 2 . 6 kilometers. The temperature soundings shows a nearly 
dry adiabatic lapse rate from the surface to an altitude of 262 meters 
and again above 805 meters. The layer between 262 meters and 805 meters 
was a stable layer which is critical to the prediction of the top of 
the surface mixing layer. 


c. Meteorological and Source Inputs to the Model 

Based on the vehicle used, in this case the Scout -Algol 
III, and the meteorological conditions at the launch time, a set of values 
for meteorological and source model input parameters for the MSFC/NASA 
Multilayer Diffusion Model were defined. Time-lapse photographs were 
used primarily as an empirical comparison between predicted and observed 
exhaust cloud rise and volume. A detailed list of the general input 
parameters required by the model is given in Appendix III. The specific 
values that were employed in the model to obtain the concentration 
mappings for the launch of the Scout-Algol III vehicle are given. 


Calculations of the maximum height (z ) of rise of the 

m 

Scout motor exhaust ground cloud, to define the location and size 
of the source, were made using rawinsonde soundings forty (40) minutes 
before launch (1430 UT) to predict the environmental dispersion of these 
exhaust effluents. The values of the ambient temperature (T) , and the 
potential temperature gradient (30 /dz) for this sounding is given in 
Table IX. This positive potential temperature gradient implies that a 
stable cloud rise relation is required (Equation A-2 or A-4). Because 
of the almost instantaneous rise of the solid fuel rockets, the instan- 
taneous rise of the solid fuel rockets, the instantaneous cloud rise 
relation (spherical entrainment) for a stable environment (Equation A-2) 
is required to describe the Scout* s ground cloud rise. The initial 
radius at the surface (r ) and the entrainment coefficient (8) in Table IX 
are required in the cloua rise calculations and are based on characteristics 
of the Scout-Algol III vehicle. Based on limited experience in predicting 
cloud rise from solid rocket motor launches at Vandenberg Air Force Base, 
this assumption appears to be justified. The time required for the cloud 
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FIGURE 12. TEMPERATURE AND DEW-POINT PROFILES FROM THE 1430 UT 

RADIOSONDE SOUNDING ON AUGUST 13, 1972, WALLOPS ISLAND 
STATION, VIRGINIA 
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FIGURE 13. WIND SPEED AND DIRECTION PROFILES FROM THE 1430 UT 
RADIOSONDE SOUNDING ON AUGUST 13, 1972, WALLOPS 
ISLAND STATION, VIRGINIA 
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TABLE IX 


THE SCOUT-ALGOL III INPUT PARAMETERS TO PLUME RISE FORMULAS 


A. General 


1. Initial cloud radius at the surface 

r 

i 

1 

0 


R 


2. Air density 

P 

167 g/m 

3. Entrainment constant 

Y I 

0.64 

4. Specific heat of air at constant 

c 

0.24 cal/g°K 

pressure 

p 



B. Meteorological Regime (1430UT, August 13, 1972) 


1. Ambient air temperature 

T 

298. 16°K 

2. Vertical potential temperature 

Be. 

0.005°K/m 

gradient 

Bz 
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where g is the gravitational acceleration, T is the ambient temperature, 
and do/dz is the potential temperature gradient. The time required for 
the ground cloud to reach stabilization height was calculated to be 482 
seconds for the parameters given in Table IX. The maximum height of 
cloud rise (Equation A-2) based on the 1430 UT sounding is 417 meters. 

This analytical prediction of the height compares well with the observed 
cloud rise heights (Figure 14) obtained using time-lapse photographs [5 ] # 
(While only a single analytically predicted cloud rise height has been 
determined, a predicted cloud rise curve can be analytically generated 
with Equation A-2.) 


The height of cloud stabilization is then employed in 
the empirical prediction of the cloud diameter based on the MSFC library 
of cloud growth observations. The expression for the dimensions .of the 
stabilized ground cloud is calculated as the standard deviation of the 
relative Lagrangian distance between two particles and is given by 


20 


yo 


= a 


xo 


v z 
T I ml 

2.15 


( 6 ) 


where a , a , and a are the vertical, crosswind and alongwind 
zo yo xo 

standard deviations for the dimensions. For the Scout-Algol III launch, 
the spherical entrainment leads directly to the above assumption of a 
symmetric stabilized ground cloud. The standard deviation for the 
predicted dimensions is 124 meters. For comparison of the predicted 
with observed cloud diameter, the photographically observed cloud growth 
as a function of time was made (Figure 15). The empirical cloud diameter, 
then, was determined to be 400 meters at an altitude of 417 meters 
(Figure 16). In the model for the source, the standard deviation for the source 
dimensions -- utilizing a Gaussian distribution -- implies that 95% of 
the material in the ground cloud was photographically visible compared to 
an expected 98%; that is, a 37o difference exists between the empirical 
results and the analytical results. 


The ground cloud became stable with the ambient air 482* 
seconds after launch and had a maximum height of about 400 meters. At 
this time, the cloud contained 1592 kilograms of the rocket engine 
effluents from approximately the first 8 seconds of engine operation. It 
will be noted that the cloud stabilization time implied by the photographic 
analysis (Figures 14 and 15) is about 2 minutes; whereas, the analytic 
results obtained using Equation 5 is about 8 minutes. The differences 
between the empirical and analytical values obtained for the stabilization 
time and dimension ar-e" explainable using the following argument. The rate 
of change in the ground clouds growth and altitude follows an exponential 
decay as the point of stabilization is approached. The empirical values 
obtained through photographic analysis are limited in accuracy; therefore, 
empricial values obtained in this manner afford only about 807o of the 
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FIGURE 16. CROSS WIND AND VERTICAL DIMENSIONS ASSIGNED 
ANALYTICALLY TO THE GROUND CLOUD FORMED FROM 
SCOUT LAUNCH AT 1430 UT, AUGUST 13, 1972, AT 
WALLOPS ISLAND, VIRGINIA 






total ground cloud development. This exponential decay of the growth 
is illustrated in the curve (Figure 17) of the rate ? s for the cloud 
radius to cloud altitude (r/z) as a function of time after launch. This 
ratio (r/z) is equal to the instantaneous entrainment coefficient (y). 

The entrainment coefficient required in the cloud rise calculations is 
the average entrainment from launch to cloud stabilization, which is the 
average value of the entrainment. The cloud stabilization can be 
defined as the time when the instantaneous entrainment coefficient 
approaches a constant limiting value. 

A detailed examination of the appropriate air temperature 
and dew point (Figure 12), and wind speed and direction (Figure 13) at 
the time of the launch of Scout-Algol III, on August 13, 1972, indicated 
the atmosphere should be treated as a single layer with regard to the 
ground cloud, rather than the nine layer model that was employed with 
Apollo 16. The NASA/MSFC Multilayer Diffusion Model 3 (spherical source) 
was applied to analyze the dispersive transport of effluents from the 
launch of the Scout-Algol III launch. (This model is summarized in the 
Section II and details for it are given in Appendix II.) 

Source and meteorological input parameters for the NASA/ 
MSFC Multilayer Diffusion Model are defined in Appendix III (Tables 
A-I and A-III) . The values of the meteorological inputs, the source 
dimensions, and the source strength in the layer were obtained relative 
to the stabilized ground cloud (Figure 16) utilizing the relations in 
Appendix III (Equation A-39 , A-40, and A-41) . These meteorological 
and source input parameters (Table X) are the basis of the effluent 
dispersion results for the launch of the Scout-Algol III vehicle on 
August 13, 1972 that were obtained with the NASA/MSFC Multilayer Diffusion 
Model. 


d. Predictions of the Centerline Concentration of Exhaust 
Effluents from the Launch of the Scout-Algol , III , 

August 13, 1973, by use of the NASA/MSFC Multilayer 
Diffusion Model 

The following predictions for the dispersion of the 
Algol III exhaust effluents from the Scout launch will focus on the 
dispersion of hydrogen chloride (HC1), carbon monoxide (CO) and alumina 
(A1 0 ) because these are the only three constituents produced in 
significant amounts that are potentially toxic and have a potential 
impact on the environmental air quality in the troposphere. While 
carbon dioxide is produced and is potentially toxic, the amount is so 
small that it can be omitted from these discussions. 

The results of the model calculations using Model 3 
of the NASA/MSFC Multilayer Diffusion Model Program and the meteorological 
and source inputs in Table X are shown in Figures 18 through 20. 
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FIGURE 17. TIME DEPENDENT ENTRAINMENT COEFFICIENT FOR THE 
SCOUT VEHICLE 


TABLE X. 


METEOROLOGICAL AND SOURCE INPUT PARAMETERS FOR THE 
SCOUT-ALGOL III LAUNCH OF AUGUST 13, 1972 


Parameter 



Model Input 


Q k - HC1 


A1 2°3 


ppm m 

r 

ppm 


m sec 


m sec 


2.153 x 10 
3. 730 x 10* 
5.033 x 10* 


RK 

a AR (T °K> 
°A TK ^ T ok) 


a = a -a 
xo yo zo 


meters 

degrees 

degrees 

degrees 

degrees 

seconds 

seconds 

meters 

dimensionless 

dimensionless 


124. 13 


meters 


meters 


degrees 

degrees 

meters 

meters 


meters 


meters 


meters 


meters 


2; 417 






CONCENTRATION (ppm) 







CO CONCENTRATION (ppm) 



FIGURE 19. MAXIMUM CENTERLINE CO CONCENTRATIONS AT THE SURFACE 

AND AT CLOUD STABILIZATION HEIGHT FOR THE SCOUT-ALGOL III 
LAUNCH OF AUGUST 13, 1972 AT WALLOPS ISLAND, VIRGINIA 
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DOWNWIND DISTANCE (kilometers) 


FIGURE 20. MAXIMUM CENTERLINE Al 2 0o CONCENTRATIONS AT THE SURFACE 

AND AT CLOUD STABILIZATION HEIGHT FOR THE SCOUT -ALGOL III 
LAUNCH OF AUGUST 13, 1972 AT WALLOPS ISLAND, VIRGINIA 






The maximum hydrogen chloride concentrations downwind 
from the Scout-Algol III launch of August 13, 1972 at a height of 
2 meters above the surface and at the cloud stabilization height of 417 
meters are given in Figure 18. Inspection of this figure shows that 
hydrogen chloride concentrations at 417 meters altitude decrease from about 2.3 
parts per million near the source to less than 0.1 parts per million at 
a distance of 6 kilometers downwind from the source. Ground-level 
hydrogen chloride concentrations reach a maximum of about 0.3 parts per 
million about 1.5 kilometers downwind from the source. At distances 
beyond 3 kilometers from the source, the hydrogen chloride concentrations 
at the cloud stabilization height and the surface are nearly equal. 

The concentrations of carbon monoxide and alumina were 
determined at the surface and at a height of 417 meters above the surface. 
Maximum concentrations of carbon monoxide at the ground surface do not 
exceed 0.52 parts per million (Figure 19) and for alumina do not exceed 
0.7 milligrams per cubic meter (Figure 20). It should be noted that the 
calculation of alumina concentrations was made under the assumption that 
the alumina emitted by the Scout-Algol III vehicle consisted of finely 
divided particulates with no appreciable settling velocities. 

The centerline concentrations at the surface and at 
the cloud stabilization altitude become approximately equal 3 kilometers 
downstream. This implies that the vertical distribution of effluents 
becomes uniform at this point and the dispersion of effluents is really 
being conducted as a Model 1 dispersion beyond this point in spite of 
the fact that Model 3 is being utilized. 

The concentration levels for this Scout launch were at 
least a factor of approximately 10 less than the public exposure limits 
(Appendix IV). Due to the small amount of effluents from the Scout, the 
isopleths of the concentrations become meaningless; therefore, they have 
not been included. 
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APPENDIX I. CLOUD RISE FORMULA [ 1] 


In order to determine whether an adiabatic or stable cloud 
rise relation should be utilized, it is necessary to determine the 
vertical potential temperature gradient / B0\ which is described by 

Vdz / 


Se. ss 9.dT + g_ = dT + g_ 

dz “ T 3z c dz c 

P P 


(A-l) 


where 0 and T are the potential and ambient temperature, g is the 
gravitational acceleration, and c^ is the specific heat of air. If 


- o 

dz 


the adiabatic cloud rise relation is used. 


(A-2) 


The maximum cloud rise z downwind from an instantaneous 
source in an adiabatic atmosphere m is given by 


ml 


2 F T 1 T 

I si 

3 2 

L 7 j u 


,41 


R 


Y- 


1/4 


R 


(A -3) 


Y t 


whereas, the maximum cloud rise z downwind from an instantaneous 
source in a stable atmosphere is given by 


ml 


8F 


LYj s 


41 


1/ 4 


+ 


R 


Y i 


R 


where F is the instantaneous buoyancy parameter / 3gQ n 


(A-4) 


, Q is the 

’ I 


4npc T. 

< p ) 

effective heat released, p is the density of ambient air, is the 

entrainment coefficient, r is the initial cloud radius at the surface, 

R 


s accounts for the vertical gradient of the potential temperature, and 

x T is the distance to reach stabilization. The subscript I means 

s i 

instantaneous and is used to flag a difference in the cloud rise models. 
The buoyancy terms, which is a function of the heat released and the 
type of entrainment, spherical and cylindrical, reflect the major 
difference in the two sources. 
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Equations (A-3) and (A-4) assume that the initial upward 
momentum imparted to the exhaust gases by reflection from the ground 
surface and launch pad hardware is insignificant in comparison with 
the thermal buoyancy flux. These relations are normally used with 
solid rocket motors. 


The following formulas for the maximum buoyant rise of clouds 
from continuous sources are also based on procedures similar to those 
given by Briggs (1970). The maximum cloud rise downwind from a 
continuous source in an adiabatic atmosphere is given by 


z 

me 


3 F x 
c sc 

* n 2 _ 3 

L 2 r c u 




(A-5) 


The maximum cloud rise z downwind from a continuous source 

me 

in a stable atmosphere is given by 


z 

me 



(A-6 ) 


where F is the continuous buoyancy flux parameter and is equal to 
c 

gQ 

c # The subscript c implies that the associated parameter is 
HpCpT 

unique to the continuous source. The primary difference in these 
continuous source relations is that the temperature constraint in the 
stable atmosphere results in a buoyancy damping. 


Equations (A-5) and (A-6) assume that the initial momentum 
flux imparted to the cloud by dynamic forces is negligible in comparison 
buoyancy flux. Again, experience in calculating cloud rise for normal 
launches of large liquid fueled rockets and for static firings has 
shown that this assumption is reasonable [6-10] . 


61 


APPENDIX II. CONCENTRATION-DOSAGE FORMULATION FOR NASA/ 
MSFC MULTILAYER DIFFUSION MODEL 


The fundamental relation for the concentration-dosage calcu- 
lation will be presented for the ellipsoidal source used in Model 3. 
These relations are appropriate to the elliptic-cylindrical distribution 
of Model 1 if the vertical dispersive interaction is neglected. This 
part of the Appendix is complex; therefore, only recommended when 
a detailed scientific knowledge is required. 


The dosage equation for Model 3 in the K-th layer is given 
by the expression 


y K’ Z BK < h<L < Z TK> = 
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r 
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L 
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r _ ( 2i ( z tk” z bk) - (^k -z k )) 2 ~j 

J 

+ exp 

2<7 zK 


+ exp 


~( 2x ( z tk~ z bk) + ( h k 2z bk +z k))' 


2 a, 


zK 


where Q corresponds to the source strength or total mass of material 
K 

in the layer and H is the height of the centroid of the stabilized 
cloud. K 


( a zK) ^ e fi- nec * ky the expression 


The standard deviation of the vertical dosage distribution 

Pk 

(A-8) 


a = a 1 x 
zK EK rzK 


X K + X zK- X rzK: (1 "% ) ' 


P K X rzK 
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where cr' describes the mean standard deviation of the wind elevation 
EK 

angle, x zK gives the vertical virtual distance, accounts for vertical 

diffusion, and x is the distance over which rectilinear vertical 
r ZK 

expansion occurs downwind from an ideal point source in the K-th layer. 

In the surface layer (K = 1), the standard deviation of the 

wind elevation angle (cr ) at the height z is described by 

ilK R 


W K=1 > 


g ER ,( Z TK^ K ~ 1 - , ) q ~( z It) q ] / T \ {A _ 9) 

(^)(^=d- R )( ZR ) q w 


where the power-law exponent (q) for the vertical profile of the 
standard deviation of the wind elevation angle in the surface layer is 


q = log I 


°ETK^ K ^ 


or. 


ER 



z tk {K=1} 


"R 


(A-10) 


here ^- s t ^ ie standard deviation of the wind elevation angle 

at the top of the surface layer. Above the surface layer (K>1), the 
standard deviation of the wind elevation angle is 


°ElJ K> ^ (°ETK + °EBK ) ( 360 ) 


(A-ll) 


where cr mTr and a are the standard deviations of the wind elevation 
ETK EBK 

angle at the top and the base of the layer. 


The vertical virtual distance x is given by the expression 

zK 


a {K} 
zo 

cr ' 

EK 


^RzK 


; a {K} s cr’ x 
* zp l J EK l'zK 


a {K} \ 1//3 K 
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X rzK\ a’ x 

' EK rzK 


-RZK + X rzK (1 ^K> •- %< K > * X rzK 


(A-12) 
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where a f K] is the standard deviation of the vertical dosage distribution 
zo L J 


at x , the distance from the source where the measurement is made in 
RzK 


the K-th layer # 


The remaining terms are common also to Model 1; that is, what 
has just been discussed is to account for the vertical expansion of 
the source cloud . 


The quantity u in Equation (A-7) is the mean cloud transport 
speed in the K-th layer. In the surface layer (K=l), the wind speed- 
height profile is defined according to the power-law expression 


°hc’ K=1 i 


{K = 1} ' P 


”r 


IC 


*R 


(A-13) 


where u is the mean wind speed measured at the reference height z and 
R R 

the power-law exponent (p) for the wind speed profile in the surface 
layer is described by 


p = log 


' 5 tk< k ” 1} 


u 


log I 


R 


' z tk< k=1 > 


‘R 


(A-14) 


here u^fK=l} corresponds to the mean wind speed at the top of the 
surface layer ^z^{K= . Thus, in the surface layer, the mean cloud 


transport speed (u[K=l} ) is 




u 


XK 


u k {K-1} - 


R 


( Z TK (&1 > - "r) V 


(z K { K =l}) 


P dz (A-15) 


"R 


which reduces to 


®E f. ( Z TK< K=1 >) UP ~( Z ll) 1 ' P - 


(A-16) 


In layers above the surface layer (K>1), the wind speed-height profile 
^u{z K ,K>lJ) is assumed linear and defined as 


u|z K , K »ij = u BK + 


' U TK U BK \ 

K z tk~ z bk) 


( z k - z bk ) 


(A-17 ) 


where u and u describe the mean wind speed at the top of the 
lis. BK 

layer and at the base of the layer, respectively. In the K-th layer 
(K>1), the mean cloud transport speed ^u^[K>l)^ is 


¥ K >i) = (« TK + a nK )/2 


(A-18) 


The standard deviation of the crosswind dosage distribution 
( cr TT i is defined by 

V yK; 


yK 


<j ' { ^ x 

AK 1 lU ryK 


+ x 


"K ‘ X yK ~ X ryK ( 1 fx) \ 
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(A-19 ) 
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+ 4.3 


*K 


-i2 
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where o' (t ) corresponds to the mean layer standard deviation of the 
AK K .1 

wind azimuth for the cloud stabilization time. ) . In the surface 

layer (K=l) , ' 

AIC { * •- > - MW&l} - ' • > 

where the standard deviation of the wind azimuth angle ^ CT * {t^ 

height z and for the cloud stabilization time t is 
R K 


at 


/ T \^/^ 

■ar^k} = <r AB { T oK} ) (ifo) 


(A-21) 
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here a A {t } is the standard deviation of the wind azimuth angle 
AR oK / v 

at height z and for the reference time period fr T \ , and the power-law 

R V oK / 

exponent (m) for the vertical profile of the standard deviation of the 
wind azimuth angle in the surface layer is 


m = log 




ATK 

cr' ( 7 k-) 

AR 1 K/ 



(A-22) 


Then 


<1tK< T K.K=1} = cr 


ATK 


{ t oK ,K=1} 



(A-23) 


where cr f t ,K=1) is the standard deviation of the wind azimuth 
ATK l oK j 

angle at the top of the surface layer for the reference time period. 
For layers above the surface (K>1), 


■W T K>K>1} - (^ TK r K } + ^bk( t k>)/ 2 < a ‘ 24 > 


where 


(£) 


(A-25) 


here cr [t } is the standard deviation of the wind azimuth angle at 
ATK ok j 

the top of the layer. 


cr' 

ABK 


r K } 



(A-26) 


here <7 } is the standard deviation of the wind azimuth angle 

ABK oK 

in degrees at the base of the layer for the reference time period 
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The crosswind virtual distance is 
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(A-27) 
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yo 


here a {K} is the standard deviation of the lateral source dimension 
yo 

in the layer at downwind distance x , x is the distance over which 

RyK ryK 

rectilinear crosswind expansion occurs downwind from an ideal point 

source, and O' describes the lateral diffusion in the layer. The 
K 


vertical wind direction shear 


K) 


in the layer is 


AO 


K 


" (° 


TK " °BK 


) (iio) 


(A-29) 


where 9 and q are the mean wind direction at the top, and at the 
TK BK 

base of the layer, respectively. 

The concentration algorithm is of the same form for the first 
three models; however, the dosage term does depend on which model 

has been utilized, and thus adjusts the concentration description to 
the specific model of interest. 

The maximum concentration for the first three models in the 
K-th layer is given by the expression 

*K {x K’ y K’ Z K> 


d k u k 

/27T <7 


xK 


(A-30) 
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where the standard deviation of the alongwind concentration distribution 
/a \ in the layer is 

V xKj 


•* ■ [('-&>)’ 


+ (J {k} 

xo . 


1/2 


(A-31) 


and the alongwind cloud length for a point source in the layer 


at the distance x from the source is 
K 


°.28(Au k )(x k ) 
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L{%} = 
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A5 K & 0 


AS K s0 


(A-32) 


where &u is the vertical wind speed shear in the layer and is defined 
K 


as 


Au k ( k=1 ) = G tk (K=1} - U R 


(A-33) 


or 


A u k (K >1} U TK “ u bk 


(A-34) 


and a {K} is the standard deviation of the alongwind source dimension 
xo 

in the layer at the point of cloud stabilization. The above equation 
for L{x } is based on the theoretical and empirical results reported 

by Tyldesley and Wallington (1967) [ 11] who analyzed ground-level 
concentration measurements made at distances of 5 to 120 kilometers 
downwind from instantaneous line-source releases. 

The maximum centerline concentration for the model in the K-th 
layer is given by the expression 


WW°’ Z K } = x K /{lateral term} (a . 35) 
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The average alongwind concentration is defined as 


X K D K /^pK 

where the ground cloud passage time in seconds is 


(A-36) 



4.3 


xK 



A-37) 


The time mean alongwind concentration in the K-th layer is defined by 
the expression 

d k 

X lJ X K ,y K ,Z K ;T A'^ = T~ 

A 

where T A is the time in seconds over which concentration is to be 
A 

averaged* The time mean alongwind concentration is equivalent to the 
average alongwind concentration when t equals T • This complex set 
of relations, then, are the computatioRs performed in Model 3 to obtain 
the concentration-dosage mappings. 
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u k t a 
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APPENDIX III. INPUT PARAMETERS FOR THE NASA/MSFC 
MULTILAYER DIFFUSION MODEL 


There are to groups of input parameters for the model. The 
source input parameters which are vehicle and meteorologically depen- 
dent (Table A-I) and the meteorological input parameters which are 
strictly dependent on meteorological conditions at launch time 
(Table A- II) . These parameters include the special set employed in 
the layer breakdown model -- Model 4. 


The source relationships given in Table A-I are determined in 


4- /v o n'k 4 14 »7 /I nl 

*»VJ> LUC DtaUJUXX^CU gJ.vui.iVk vxvuU. 


T’Vic 


c* f* a r\A a irr\ Ac 


sxi* A t 


r\ f 


the crosswind source is 
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= y£k! 
4.3 


(A-39) 


and the standard deviation of the alongwind source is 


a 



= 

4.3 


The source strength in the K-th layer is 


(A-40) 
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YfK] Xf KV 
S Y{K} X{K} 


Q 

T 


z -z 
TK BK 


(A-41) 


where Y{K} and X{K} describe the crosswind and alongwind dimensions of 
the cloud in the K-th layer, and Q is the total source strength in the 

ground cloud in units of mass. 


Equations (A-39) and (A-40) are based on the assumption that 
the alongwind and crosswind distribution of material in each layer is 
Gaussian and that the visible edge of the cloud represents the point 
at which the concentration is one-tenth the concentration at the cloud 
center in the K-th layer. Equation (A-41) assumes that the cloud is 
spheroidal in the plane of the horizon and that the total source 
strength in the K-th layer is given by the relative cloud volume in 
the K-th layer. Because the models require the source strength per 
unit height, the total source strength in the K-th layer must be divided 
by the depth of the layer. 

The first nine meteorological parameters follow directly from 
the thermodynamic and kinematic profiles of the atmosphere. The remain- 
ing two parameters (layer model) are empirical atmospheric constants. 
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TABLE A-I 


SOURCE INPUTS FOR THE MULTILAYER MODEL CALCULATIONS 


PARAMETER 


Layer 
Model 1, 
2, 3 

Layer Break- 
down Model: 4 

DEFINITION 

z 

R 

z 

R 

Reference height in the surface layer 

z 

z 

Height of the layer base 

BK 

BL 

z 

TK 


Height of the layer top 

T 

T 

Source (cloud) stabilization time 

K 

L 


X 

X 

Distance over which rectilinear lateral 

ryK 

ryL 

expansion occurs downwind from an ideal 
point source 

a {K5 


Standard deviation of the crosswind source 

yo 


dimension in the K-th layer 

a {K} 


Standard deviation of the alongwind source 

xo 


dimension in the K-th layer 


t* 

Time of layer breakdown 

Q 

K 


Source strength in the layer 

J 


Scaling coefficient 
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TABLE A- II 


LIST OF METEOROLOGICAL MODEL INPUTS 


PARAMETER 


Layer 
Model: 1, 
2, 3 

Layer Break- 
down Model: 4 

DEFINITION 

^R 

u 

RL 

Mean wind speed at reference height z 

R 

u 

u 

Mean wind speed at the base of the layer 

BK 

BL 


u 

U 

Mean wind speed at the top of the layer 

TK 

TL 

e 

e 

Mean wind direction at the base of the layer 

BK 

BL 


e 

9 

Mean wind direction at the top of the layer 

TK 

TL 

a {t } 

a [t } 

Standard deviation of the wind azimuth angle 

AR or 

ARL OL 

at height z 

R 

a ABK^ T oK^ 

CT ABL^ T oL^ 

Standard deviation of the wind azimuth angle 

at the base of the layer 

a {t } 
ATK or 

a {t } 

Standard deviation of the wind azimuth angle 

ATL oL 

at the top of the layer 

T 

T 

Reference time period 

oK 

oL 

a- 

O' 

Lateral diffusion coefficient 

K 

L 


p{K=l} 

p{L»l} 

Power-law exponent of the wind speed profile 

a 

in the surface layer 

Standard deviation of the wind elevation 


ERL 

angle at height z 

R 


o__ 

Standard deviation of the wind elevation 


EBL 

angle at the base of the L-th layer 


a 

Standard deviation of the wind elevation 


ETL 

angle at the top of the L-th layer 



Vertical diffusion coefficient 
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APPENDIX IV. TOXICITY CRITERIA 


Realistic evaluation of the potential 'hazard arising from 
high near-field concentrations of toxic effluents from solid rocket 
exhaust requires both a knowledge of the surface deposition of these 
effluents -- which can be obtained with the MSFC/NASA Multilayer 
Diffusion Model (Appendix II), and a toxicity criteria to evaluate the 
hazard from this surface deposition of effluent — which is the incum- 
bency for this discussion. The Federal Air Quality Criteria does not 
presently include any of the liquid or solid rocket exhaust effluents; 
however, the National Academy of Sciences does afford definite guide- 
lines for the exposure to the toxic effluents associated with thdse 
exhausts. These guidelines are ecologically sound, based on the current 
limited knowledge of the effects of these effluents, and are the basis 
of the toxicity criteria that will be* given [12, 13]. 

The primary effluents from any solid rocket exhaust are: 
aluminum oxide, (A1 0^), hydrogen chloride (HC1) , carbon monoxide (CO), 
carbon dioxide (CO 2 ; , hydrogen (H 2 ), nitrogen (^) and water vapor 
(H 2 O). While only the first four compounds are toxic in significant 
concentrations, there is always a potential hazard of suffocation 
from any gas which results in the reduction of the partial pressure of 
oxygen to a level below 135 mm Hg (18% by volume at STP) . Oxygen level 
reduction does not appear to be a hazard from rocket exhaust due to the 
large volume of air which is entrained into these exhaust clouds; 
therefore, this potential hazard can be neglected in this discussion and 
the attention directed to only the initial four toxic compounds. (A 
liquid rocket motor has only one toxic effluent — carbon monoxide.) 

The exposure level for toxic effluents are divided into three 
categories: public exposure level, emergency public exposure level, 
and occupational exposure level. The public exposure levels are 
designed to prevent any detrimental health effects both to all classes 
of human beings (children, men, women, the elderly, those of poor 
healt, etc.) and to all forms of biological life. The emergency level 
is designed as a limit in which some detrimental effects may occur, 
especially, to biological life. The occupational level gives the 
maximum allowable concentration which a man in good health can tolerate 
this level could be hazardous to various forms of biological life. 

The toxicity criteria for the toxic effluents in solid rocket 
exhausts are given in Table A-III. Public health levels for aluminum 
oxide are not given because the experience with these particulate is 
so limited that, at best, the industrial limits are just good estimates. 
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TABLE A-III. AIR, QUALITY TOXICITY STANDARDS** 


TOXIC SOLID 
ROCKET EXHAUST 
PRODUCT 

TIME 

INTERVAL 

(MINUTES) 

CONCENTRATION 

PUBLIC 

EMERGENCY 

OCCUPATIONAL 



*** 





5.0 mg/m;? 

X 

50 mg/m;? 

Alumina (Al n 0_) 

30 

2.5 mg /nr 

X 

25 mg/m' 5 

2 3 

60 

1.5 mg /nr 

X 

15 mg /nr? 


480 

1.0 mg/nr 

X 

10 mg/nr 


10 

wEsm 

BB 

30 ppm 

Hydrogen Chloride 

30 


mmSm 

20 ppm 

(HC1) 

60 

BB 

H9 

10 ppm 


10 

90 ppm 

i 

275 ppm 

1000 (1500*) ppm 

Carbon Monoxide 

30 

35 ppm 

100 ppm 

500 (800*) ppm 

(CO) 

1 

60 

25 ppm 

66 ppm 

200 (400*) ppm 

DOSAGE: 


200 ppm/ 





time interval 



Carbon Dioxide 

480 

X 

X 

Average - 5000 ppm 

(co 2 ) 


X 

X 

Peak - 625 0 ppm 


Parts of vapor or gas per million parts of contaminated air by volume at 
25°C and 760 mm Hg. 


**These values were reviewed on the phone by Ralph C. Wands, Director 
Advisory Center on Toxicology, National Academy of Sciences, 2101 
Constitution Avenue, Washington, D. C. 20418, April 1973. 

*At these concentrations, headaches will occur along with a loss in 
work efficiency. 

*** EPA suggests a safety factor of ten (10) to be applied to occupational 
exposure limits [ 
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Hydrogen chloride is an irritant; therefore, the concentration 
criterion for an interval should not be exceeded [13], Since hydrogen 
chloride is detrimental to biological life, and in view of the fact 
that most launch sites are encompassed by wild life refuges, the 
emergency and industrial criteria for hydrogen chloride are not appro- 
priate to the ecological constraints. Because of the large volume of 
air entrained in the exhaust cloud, the potential hazard from carbon 
monoxide and carbon dioxide can be, in general, neglected. 

Any detrimental health effects due to combined toxicological 
action of these ingredients has been omitted because of a lack of 
knowledge in this area. However, investigations are currently underway 
to study this problem and to learn more about the biological effects of 
hydrogen chloride. 
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